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In an analysis of the US, the UK, and the German stock market we find a change in the behavior 
based on the stock’s beta values. Before 2006 risky trades were concentrated on stocks in the IT 
and technology sector. Afterwards risky trading takes place for stocks from the financial sector. 
We show that an agent-based model can reproduce these changes. We further show that the 
initial impulse for the transition might stem from the increase of high frequency trading at that time. 
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I. INTRODUCTION 

In this paper we analyze transitions in the structure 
of the US, the UK, and the German stock market. In 
particular we observe a phase of dominance of IT oriented 
stocks followed by a transition period that leads to a 
dominance of the financial sector. 

The analysis of the differences in the returns of stocks 
have long been dominated by the discussions around dif¬ 
ferent versions of a CAPM model [l^, [2^. The original 
version of the CAPM is in fact a one factor model, which 
postulates that the returns of the stocks should be 
governed by the market return tm and only differ by an 
idiosyncratic component Pi for each stock z, such that 

rpt) = PiTmit) + u(t). (1) 

Hence, stocks differ by the amount of volatility with re¬ 
spect to the market, and economic rational necessitates 
that higher stock volatility is compensated by higher ab¬ 
solute returns (additionally eq. may incorporate the risk 
free interest rate). Empirical tests of this model had 
rather mixed results and have let to two conclusions: 
More factors are needed to explain the variation of stock 
returns. The widely used Fama and French [T^ model 
for example is a three-factor model that incorporates firm 
size and book-to-market ratio. The second conclusion 
was that beta values are not constant but time-varying, 
see Q- The reasons for the variability of the betas are 
manifold. They could change due to microeconomic fac¬ 
tors, the business environment, macroeconomic factors, 
or due to changes of expectations, see, e.g., Qilll- Also 
models that assume a first-order auto-regressive process 
have been suggested, see 

Our approach to identify different states of a stock 
market consists in an analysis of a covariance matrix, 
similar to [l^, and of the transaction volumes, like in 


. The properties of the covariance matrix of asset re¬ 
turns depend on the time horizon T in which they are 
determined. For short T in the order of months they are 
rather volatile, and partly mirror economic and political 
changes [H, [ig . Q for example argue that correlations 
increase in times of crisis, which has profound implication 
for portfolio choice and hedging of risks. 

For large T in the order of several years, a principal 
component analysis [l3, of the correlation matrix 

is possible. The component with the largest eigenvalue 
can be interpreted as the market. The P coefficients are 
proportional to the corresponding eigenvector. 

In order to detect changes in stock betas we use time 
windows of less than 4 years and a rather large numbers 
N of stocks for different markets. In this case the prin¬ 
cipal component is well separated from the rest. Within 
the assumption of market dominance motivated by imple¬ 
menting eq. with a stochastic volatility model (SVM), 
one can determine the P coefficients. A problem may be 
the statistical accuracy, which could be of order \JN/T 
as suggested by random matrix theory (^ . However, a 
Monte Carlo simulat ion s hows that the errors for the P 
are in the order of \p2PT. 

This paper consists in a substantial extension of the 
research presented in [2^ , with respect both to methods 
and data. The paper is organized as follows: In section 
m we briefly describe the data sets before we describe 
the methodology to analyze the covariance matrix and 
the distribution of the stock returns. In section m we 
present a Monte Carlo study for the error estimate on 
p. After this we show the transitions in the markets and 
introduce a sector specific risk measure. In section El we 
present a model that can replicate the transitions and we 
discuss whether the cause of the changes is an internal or 
external one. Section EH concludes. 
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US 

S&P500 

UK 

FTSE350 

Germany 

CDAX 

period 

1995-2013 

1997-2013 

1999-2013 

T 

4782 

4294 

3691 

N 

356 

132 

78 

sector 




Energy 

32 

5 

0 

Materials 

23 

7 

7 

Industrials 

51 

30 

25 

Cons. Discr. 

56 

26 

12 

Cons. Staples 

35 

10 

5 

Health 

32 

5 

10 

Financial 

60 

37 

7 

Technology/IT 

37 

4 

7 

Telecom. 

3 

4 

3 

Utilities 

27 

4 

2 


TABLE I. Summary statistics of the data sets 


II. MATERIALS AND METHODS 

A. The data sets 

For our analysis we use data from Thompson Reuters 
on the closing price of stocks which were continuously 
traded with sufficient volume throughout the sample pe¬ 
riod and had a meaningful market capitalization [? ]. For 
the US we choose stocks which are part of the S&P500 
stock index. For the UK the stocks in our sample are 
listed in the FTSE350, the German stocks are all part of 
the CD AX (and are with very few exemptions also listed 
in the MDAX, SDAX or DAX30). The size of the US 
market allows us to collect a time series corresponding 
to 20 years of data. For the European markets it is not 
possible to analyze a quite as long time horizon, since not 
enough stocks have been traded for such a time span. We 
have collected the sector classification of the firms, using 
the GIGS classification for the US market and the (for our 
purposes practically identical) TRBC classification from 
Thompson Reuters for the European markets. Table H] 
summarizes the data sets and the sector information. 


B. Analysis Method for Correlations 

Stock markets can be analyzed by the study of the 
correlation between the returns of the participating 
firms. The N firms are indexed by i = I,-- - ,N. 
A return is given by the log of the price ratio be¬ 
tween consecutive days. The returns are normalized by 
J2t=1,To Ei=i.iv rfir) = NTq. t denotes the days in one 
time window. For the covariance matrix C we consider 
time windows of size T centered at time t. C is given by 


with the abbreviation for the time average 

T<t+Tl2 

= E (3) 

T=t-T/2 

for any observable A. In eq. © the small time averages 
< Ti > are neglected. When C is derived from the returns 
of many stocks in a long time window T oc Tq, one usually 
observes that the matrix C has one large eigenvalue Aq in 
the order of N with a corresponding eigenvector that we 
denote /3i. All Pi have the same sign and can be chosen 
positive. We normalize by Pf = N. The remaining 
eigenvalues are of order of 1. The first eigenvector can, for 
example within the framework of a principal component 
analysis, be interpreted as the market. This means that 
this eigenvector can be interpreted as the weights of the 
single stocks within the market factor. Hence, a market 
return tm can be defined by the the projection of r on /3 

rnir) = ^Y^PirPr) (4) 


Due to the relation 


{nrM)To,t 


the components of the leading eigenvector are P- 
coefhcients in a CAPM approach (leaving out the risk¬ 
free interest rate). With T = Tq we would have only one 
vector Pi centered at time {To/2). A time dependence of 
P can be achieved by using a moderate time window T 
(in the order of years). 

To derive meaningful Ps we assume that the return 
follows a stochastic volatility model (see, e.g., [3, H^): 
The returns are the product of a noise factor and a slowly 
varying stochastic volatility factor. The latter should be 
considered as constant over the window size T. Then 
eq. ([5]) corresponds to an average over the noise with a 
statistical error depending on the properties of r^. 

As a first example we consider rpr) = 'yrjiT with an 
i.i.d. Gaussian noise rj. For a finite T we obtain a 
Marcenko-Pastur spectrum spread over an interval 
7^(1 ± 2yjN/T) (instead of the degenerate eigenvalue 
7 ^ ). For N ^ 400 a time window of only a few years 
would lead to prohibitive large uncertainty. However, 
this model cannot account for the occurrence of one large 
eigenvalue. 

This can be reproduced by the second example with 
ri{T) = 'Yi'rjT- In this model all stocks follow the market 
described by Gaussian noise. For T —^ oo the covariance 
matrix C has one eigenvalue Ao = if with eigenvector 
Pi oc 7 i and N — 1 zero eigenvalues. At finite T the 
eigenvectors and the zero eigenvalues are unchanged. Aq 
is multiplied with a distributed number with mean 1 
and variance 2/T. To describe the observed spectrum of 
small eigenvalues we consider a second process that leads 
to an additional additive component Cuj in C. 
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We assume market dominance in the sense that 7 ^ is of 
order N and ( 7 , with constants is of or¬ 

der 1. Perturbation theory for large N, see the appendix, 
shows that Ci does not change Aq and /3i up to 1/N 
contributions. The remaining eigenvalues are strongly 
dependent on the noise. Only their sum is given by 
trace(C)-Ao. Neglecting very small quantities, < >T,t 
is a measure of the volatility in the window. 

Ao determines the size of the market return < > via 

AA, = ^ ( 7 ) 

C. The shape parameter of the returns distribution 

In order to analyze changes in the distribution of the 
stock returns we estimate the tail parameter of its pdf 

S ). We characterize / by a Pareto-Feller distribution 
, where / depends only on and a hnite /(O). The 
two parameters are a scale parameter ro and a tail index 
a. It is given by 



FIG. 1. log(l-|-N(A)) for the empirical spectrum (histogram) 
N{\) from S&P in 2004 and simulated spectrum (line). For 
the simulation we used log-normal distributed j3i with mean 
0.93, 9 = 0.26 and log normal distributed 7 ii with 70 = 0.865. 


N 

log-norm. 

Laplace normal ju — 1 

T 

log-norm. 

100 

0.109 

0.109 

0.108 

0.125 

500 

0.131 

200 

0.107 

0.109 

0.110 

0.124 

1000 

0.094 

400 

0.106 

0.108 

0.110 

0.123 

2000 

0.066 


/(r) oc 



{a - 2)rl 


-{a+l)/2 


( 8 ) 


Performing fits with limited statistics a and tq are 
strongly correlated. Therefore we fix ro by the condi¬ 
tion ro = E[r^]. 


III. MONTE CARLO SIMULATIONS 

From the discussion in section fll Bl we expect an error 
on the leading eigenvalue Ao in the order of \/2/T. This 
does not imply the same accuracy for the eigenvector /3. 
To estimate the size of errors in windows of several years, 
we perform a Monte Carlo study based on a fairly general 
SVM. In a single window we assume the following returns 

ri{t) = Ve 7 oi r]t + Vl-O -fu Vit (9) 

Comparing equation ® with the CAPM definition (jH we 
see that the first term corresponds to the market com¬ 
ponent with strength 9 and the second term can be in¬ 
terpreted as the idiosyncratic component due to trading 
activity for specific stocks. For the i.i.d. noise factors we 
use r]t iV(0,1) and rja N{0, 1). We checked that a Lapla- 
cian noise for rju as suggested by 0 does not change the 
result. The parameters 7 oi and ju are independent of 
time and normalized to 7 ^^ = N. For a window size 
T —>• 00 we get for C 


( 10 ) 


TABLE 11. The left part of the table shows the average error 
(5% confidence level) of /3 at T = 750 for various N and 
models for 71 . The right shows the error for N = 356 and 
log-normally distributed 71 for different T. 


Using perturbation theory (see the appendix) the leading 
eigenvalue and its eigenvector jSi = y/Nfi are given up 
to terms of order of 1/N by 

Ao = 9{N - 1) -h 1 and A = 70 * (11) 

We simulate the returns from eq. m with given values 
for 9, Pi and 71 ^ in a finite window. From the eigenvectors 
of the simulated covariance matrix we can estimate the 
statistical error on P due to the finite T. For this we 
need reliable values of the input parameters. The market 
strength 9 follows from the well measured empirical Aq. 
For the input Pi we use a log-normal distribution, which 
represents the observed spectrum very well. For 'ju we 
choose a normal, Laplace and log-normal distribution, 
which depend on two parameters. Since the mean of 7 ^ 
must be 1 only the mean 70 is a free parameter. 

We determine 70 by optimizing the agreement of the 
observed eigenvalue spectrum N(X) with the simulated 
spectrum. As an example we show this comparison us¬ 
ing log{l+N{X)) for the S&P data in a window of 3 years 
around 2004 and log-normally distributed 71 . Both agree 
surprisingly well. A Kolmogorov-Smirnov test leads to 
a p-value of 0.2. In contrast, Laplace or normally dis¬ 
tributed 7 ii lead to p ^ 10“^ or less. The model ac¬ 
counts apart from the bulk also for the isolated medium 


Cij = 6»7o*7oj + (I - 6»)7ij (5y 
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rank(|3) 


FIG. 2. Input p as function of the rank together with the 
simulated 5% confidence level error range for various window 
sizes T and N — 356. 71 are log-normal distributed with 
70 = 0.865 and 6 = 0.26. Plotted with an offset of 1 for each 
series. 


eigenvalues attributed in the literature to sub-markets 
. Only the second largest eigenvalue attributed to 
the trading volume [ 2 ^ does not correspond to a statis¬ 
tical fluctuation. 

The result of 200 Monte Carlo repetitions of the dy¬ 
namic eq. © is shown in figure [5] for various T. The 
ordered input j3i are connected by a line. The errors 
correspond to a 5% confidential range for a single mea¬ 
surement. There is little dependence on N or the as¬ 
sumed distribution of 71 . In table HIl we give the average 
error for various N and different 71 and T. The errors 
only vary with T by 1/ \/T. These simulations prove that 
within the assumed SVM the values of Pi can be reliably 
estimated also for moderate window sizes T. 

To summarize, for the empirical analysis of C in the 
next section we make the following assumptions: From 
the market hypothesis we can establish the leading eigen¬ 
vector of C as CAPM /3-coefficients. By the SVM as¬ 
sumption the time average in eq. m corresponds to an 
average over the noise. Making the market dominance 
assumption the errors on Aq and Pi are of the order of 
1/N, v^. 


IV. TRANSITION OF THE MARKETS IN 2006 

We apply our approach to 356 stocks from the S&P 
market, 132 stocks of the British FTSE market, and to 78 
stocks from the German market. To obtain the possible 
minimum window size T we look at the large eigenvalue 
Ao(t) and the corresponding eigenvector. As the criterion 
we use the presence of (only) positive values of Pi{t). In 
this way we find for T a value of roughly 3-4 years for all 
markets. For a better visualization of the time variation 
we use overlapping windows by varying t in steps of years. 


3.5 
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1996 1998 2000 2002 2004 2006 2008 2010 2012 2014 

FIG. 3. Time dependence of Pi for 356 stocks of the S&P 
market. The 35 stocks with largest P in 1998-2002 are shown 
in red, the 20 largest in 2007-2010 in blue. 



FIG. 4. Time dependence of Pi of the British FTSE market. 
The 7 stocks with largest P in 1998-2002 are shown in red, 
the 7 largest in 2007-2010 in blue. 

In figure[3]we show the /3—coefficients derived from the 
largest eigenvector of C for the S&P market for a time 
window of 3 years. Except one case around 2001 they 
are all positive. Some of the stocks exhibit a substantial 
time variation with a transition around 2006. Stocks with 
large P during the years 1998-2002 (this time interval is 
called ITB for IT bubble hereafter) change to small P 
values around 2006, their values remain low in 2007-2010 
(this time interval is called FB for the finance bubble 
hereafter). Vice versa those stocks with a large P in the 
finance bubble exhibit small values before 2006. A similar 
effect occurs also for the FTSE market (shown in figure 
|4]) and the German market (shown in figure [5]). For both 
a window size of 4 years is used. 

R{t, s)=AsY, - 1.0)ft (t) V{t, i) (12) 

ies 

The normalization constant As is chosen to have 
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FIG. 5. Time dependence of j3i for 78 stocks of the German 
DAX market.The 15 stocks with largest /3 in 1998-2002 are 
shown in red, the 15 largest in 2007-2010 in blue. 



FIG. 6. Time dependence of the risk parameter R(t, s) for 
the eight sectors with R ^ 0 of the S&P market. 


A more detailed characterization of the market can 
be obtained by considering the sector s out of the 
GICS/TRBS classification for all firms. An inspection of 
the firms with large /3 during the ITB in figure [3] shows 
that they dominantly belong to the IT/technology sec¬ 
tor. Likewise firms with large /3 during the FB are mostly 
from the financial sector. Since a ,5 > 1 signals a risky in¬ 
vestment, we can define a market risk measure R{t, s) for 
the sectors by multiplying Pi > 1 with the number V{t,i) 
of traded shares in each window. Note that for the fol¬ 
lowing analysis we merge the sectors IT and telecommu¬ 
nication for the UK and Germany since we have only few 
stocks in these sectors and they show similar behavior. 

In figure [S] the risk parameters from eq. (1121) for the 
S&P market is shown as a function of time. Only the 
technology sector (red) before the transition in 2006 and 
the financial sector (blue) after 2006 exhibit large values 
of the risk measure. The value of the risk measure is 



FIG. 7. Time dependence of the risk parameter R{t, s) for 
the nine sectors of the FTSE market. 



FIG. 8. Time dependence of the risk parameter R{t, s) for 
the nine sectors of the DAX market. 


small for all other sectors. Due to the time window of 
3 years the time of the transition can be fixed only with 
an error of 1.5 years. A similar phenomenon is seen for 
the FTSE market in figure Hand the German market in 
figure [51 In contrast to the other markets the industrial 
sector of the DAX shows a peak in R at 2006. This may 
be due to the fact that this sector contains one third 
of all firms. Some are large firms that are difficult to 
pinpoint to a specific sector. Since R is normalized to 
J2s = 1) the small R of technology and financials 

are compensated mainly by the industrial sector. The 
transition for the S&P appears to be somewhat sharper 
than for FTSE and DAX due to the smaller number of 
stocks in the latter. 
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FIG. 9. (a) Time dependence of response time and traded volume at the Eurex exchange, source: Eurex Exchange [H. (b) 
Normalized monthly volatility of the S&P, FTSE and DAX indices, normalized VIX index. 



FIG. 10. Simulation of R as function of time with the diluted 
Ising model. Parameters are g = 3.99, A — 100 and S = 8. 
For 100 < t < 120 a field h — —0.03 is applied. The black 
line gives the average of i?(0). 


V. EXPLANATIONS FOR THE OBSERVED 
TRANSITION 

A. Interpretation with a diluted Ising model 

The behaviour of the risk parameter R in figures |B][5] 
indicates a phase transition analogous to models in sta¬ 
tistical physics with R as order parameter. Such a model 
can be constructed by a generalization of the Ising model. 
There are A agents trading one stock per time out of the 
S sectors. Each agent a is characterized by a spin value 
ctq. Values Ca = ±1 denote trading in the risky sectors 
IT or financials, and (Jq = 0 denotes the remaining S — 2 
sectors. We assume that the /3-dependent factor in eq. 
m can be replaced by its mean. Then the normalized 
R{s) is equal to the fraction of agents trading in sector 


s. The agents can change their opinion due to an inter¬ 
action between all other agents. At each time they chose 
a new value of (Tq by the following probabilities w{a) 


i(;(±l) = — exp{±g{m + h)) 

(13) 

Wn 


,^(0) = ^(5--2) 

(14) 


where m A = 'Yha'^o. and g is the strength of the in¬ 
teraction. Wn normalizes the probabilities, h denotes a 
possible external field. Since agents with (Ta = 0 do not 
contribute to m, the model corresponds to a dynamical 
dilution. For small g the system is in the disordered state 
with i?(s) = 1/5' and for large g in the ordered state with 
one of the i?(±l) becoming large. 

The model can be solved analytically for large agent 
numbers A, as shown in appendix m For 5 > 6 a first 
order transition occurs at a critical value gc- An internal 
reason for the transition can be modeled by a dynami¬ 
cally changing g. However, this faces the following prob¬ 
lem: To observe a constant R over 6-10 years, like for the 
S&P in figure [6l the time constant to change g must be 
in the same order of magnitude. In the two years around 
2006 the S&P market changes from an ordered state into 
a disordered and back into another ordered state, making 
an external reason more likely. 

Reviewing possible external events around 2006 there 
seems to be no major political event nor any drastic 
change in asset prices. In fact, the volatility around that 
time was relatively low, as shown in the right panel of 
figure [9l One observes similarly high volatility before 
and after the transition. The only event seems to be the 
onset of high frequency trading (HFT) in 2005, see, e.g., 
[iEini- In the left panel of figure[5]we show (as a rep¬ 
resentative example) the response time and the traded 
volume at the Eurex exchange. The small response time 
and a maximum of the trading volume hint at a growing 
dominance of computerized HFT trading after 2005. 
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FIG. 11. Time dependence of the shape parameter a. For bet¬ 
ter readability we added 3 (5) to the values of DAX (FTSE). 

There are several possibilities how HFT can trigger 
the transition. Faster trading may reduce the time con¬ 
stant in the dynamics of g. Generally, computer pro¬ 
grams might be more able than humans in hedging risks. 
Before 2008 they used almost riskless strategies, as arbi¬ 
trage or flash trading, which became less important af¬ 
terwards. In the model these effects can be accounted for 
by a change of g or the effect of a field h. If one chooses 
g near g^, both phases are coexistent and only very small 
changes of g or are needed to change the phase. In fig¬ 
ure [ini we show a simulation of R with the probabilities 
from eq. 0 with constant g and application of a small 
field h at 100 < t < 120, which disfavors a previous risky 
sector. Obviously the model can reproduce the observed 
R for the markets. 


B. High frequency trading and the returns 
distribution 

The appearance of HFT should leave traces in the dis¬ 
tribution of returns. Advocates of HFT m claim that it 
leads to a more efficient market. More efficiency should 
lead to less price changes and therefore to an excess of 
smaller returns. Critics [2^ assert that computerized 
trading increases instabilities, which amounts to larger 
returns. Both effects can be seen in the pdf for the market 
return tm ■ Its pdf can be characterized by the shape pa¬ 
rameter a using the Pareto-Feller parametrization from 
eq. We obtain a by maximizing the Log-Likelihood 
L in each window. Errors on a correspond to a change 
of L by 0.5. 

In figure [TT] we show the time dependence of a for the 
three markets. Before 2006 one finds values a ^ 4 — 5 
with good probabilities. For all three markets a drop 
to values below 3 appears after 2006. Small a imply a 
much more enhanced tail of the pdf as expected from 


HFT. The x^ probabilities are worse in 2006, but still 
acceptable on the 5% level. However, the lower probabil¬ 
ities are due to systematic deviations from ([8|) . 



FIG. 12. Distribution / of the standardized (< r\j >= 1) 
market return \rM\ for the S&P market. For \rM > 1.5| / is 
multiplied with 20. 

In figure [12] we show some typical pdfs of the market 
returns before and after 2006 for the S&P market. We 
see a perfect description by eq. for the returns from 
the time windows centered in 2000 and 2002, whereas 
for those centered in 2008 and 2009 a substantial excess 
at r ^ 0 occurs, and the badly described tail extend to 
much larger value as before. This behavior is expected 
from HFT. 

For the European markets, shown in figures |T3| and 
M in the appendix, the excess of small returns is less 
significant. Only a significantly enlarged tail is observed 
after 2006. These markets might be less affected by HFT 
and therefore only the instability effect is seen. There 
is no trace in the volatility, since both effects can cancel 
out in E[r^]. 


VI. CONCLUSIONS 

The literature on regularities in asset returns has for 
a long time argued that the f3 values of stocks are time 
varying. We have shown that we can extend the con¬ 
cept of j3 values to systematically describe a risk measure 
for stocks from different sectors of the economy. This 
slowly varying sector specific risk measure describes or¬ 
dered states in the market and identifies sectors which 
show concentration of market risk. A possible trigger 
for the observed transition may be the onset of high fre¬ 
quency trading in 2005. 
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Appendix A: Perturbation Theory 


Assume a matrix C can be written as C = Cq+Ci with 
a small perturbation Ci. Cq has one large eigenvalue 
Eq and N — 1 degenerate zero eigenvalues. Due to the 
degeneracy we can impose for the eigenvectors with 
/i > 0 of Co the conditions 


e'' ■ (Cl e^) = 0 for /r, i/ > 0, i/ yf /i (Al) 


The eigenvalues and eigenvectors /f of C can be ex¬ 
panded in a power serie in Ci/Aq. For ^ = 0 we get 


Ao — Eq + c^-{C l e°)J- — 
Eq 


e°-(C^°)-(e°.(Cie°))^' 

(A2) 


/r = 




(Cie° 


Eo 


(Cie° 


(A3) 


The remaining eigenvectors need the solution of condition 

(ED 


A^ = e'^.(Cie^)-^(e°.(Cie'^))^ (A4) 

/f = er--^(e°-(Cie'^)ef (A5) 

Eq 

Note that this expansion reproduces the exact result for 
tr C and for Ci proportional to a unit matrix. 
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With (C'o)y = 7 i 7 j we have Eq = ( 7 , 7 ) = 7 ^ and e° = 
Inserting ( 7 ,( 717 ) = Ak 7 ^ into eqns. (IA2I) and 
(IA3I) we get for Aq and Pi = \/]V/° 

Aq = 7^ + Ai H—2(^2 ~ ^1) (^6) 

7 



with = (iV/ 7 ^)A 2 . Market dominance implies i?o = 
7 ^ oc -/V and the constants Aj; are of order 1. Eqns. 
(IMl) and (ITfll show that up to corrections of order 1/N 
the leading eigenvalue Aq and its eigenvector Pi do not 
depend on Ci. 


Appendix B: Diluted Ising Model 


The transition probabilities (US correspond to the 
heat bath algorithm for the following equilibrium distri¬ 
bution 

w{a) = — exp [a • g{ra(a'f‘ 12 -|- /iTO(cr))] (Bl) 

Zj 


with m(cr) = (1/A) The partitioned sum Z 

we calculate by using the Gaus trick y/n eyq){mp) = 
f dx exp(—a;^-|- 2 ma;) and evaluating the integral for large 
A. We get for Z 


InZ = A 


ln(S' — 2 + 2ch{{mo + h)g)) — ^TOq 


(B2) 


The expectation value uiq of m{a) must maximize InZ. 
This leads to the so called mean field condition for the 
order parameter mg 


Too 


2sh{mQ{g + h)) 

S — 2 + 2ch{mo{g + h)) 


(B3) 


The fraction at = 0 of agents with a a = ±1 is given 
by 


R± 


exp(±mog) 

S' — 2 -I- 2ch{mog) 


(B4) 


Eqns. (IB^ and (IbHI) lead already at = 0 to a surpris¬ 
ingly rich spectrum of phases depending on the number 
of sectors S and g. For S < 6 we have a similar behaviour 
as in the Ising model (S = 2). At gc = S/2 a second or¬ 
der transition occurs. At S = 6 the transition is still of 
second order, but with different critical exponents. The 
general case S > 6 is illustrated by a plot of — \nZ/[gA) 
in figure 1131 Below g < gi only the disordered phase ex¬ 
ists. At g = gi a toq 7 ^ 0 solution appears corresponding 
a metastable ordered phase, since InZ < ln.Z(mo = 0). 
At g = gc both phases coexist. For g > gc the mg = 0 
phase becomes metastable. Finally for 5 > <?2 = S/2 
only the mg 7 ^ 0 solution exists. The critical values of g 
and the values of R{1) at criticality corresponding to a 
lower bound are given in table IIIII 



FIG. 13. — In Z/{gA) as function of a; = gmo for various 

values of g. Values for tj < 1 are multiplied with 10. 


s 

gi 

5c 

52 

Rig.) 

6 

- 

3 

- 

1/6 

8 

3.73 

3.82 

4 

0.70 

9 

3.97 

4.20 

9/2 

0.81 

10 

4.19 

4.58 

5 

0.87 


TABLE III. Critical values of g and R{gc) at g = gc. 


Appendix C: Return distributions for the European 
markets 



|r| 


FIG. 14. Distribution of / for the standardized (< r%f >= 1) 
market return \rM\ for the FTSE market. For |rM > 1.1| / is 
multiplied with 10. 
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FIG. 15. Distribution of / for the standardized (< r\f >— 1) 
market return |rM| for the DAX market. For \rM > 1.1| / is 
multiplied with 10. 







